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Abstract. Measured hard photon multiplicities from second-chance nucleon-nucleon collisions are used 
in combination with a kinetic thermal model, to estimate the break-up times of excited nuclear systems 
produced in nucleus-nucleus reactions at intermediate energies. The obtained nuclear break-up time for the 
129 Xe + nat Sn reaction at 50,4 MeV is At « 100 - 300 fm/c for all reaction centralities. The lifetime of the 
radiating sources produced in seven other different heavy-ion reactions studied by the TAPS experiment 
are consistent with At ~ 100 fm/c, such relatively long thermal photon emission times do not support the 
interpretation of nuclear breakup as due to a fast spinodal process for the heavy nuclear systems studied. 

PACS. 21.65.+f - 25.70.-z - 13.75.Cs - 64.70.Dv 



>"H , 1 Introduction 
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Nucleus-nucleus collisions are the only available tool to 
explore in the laboratory the different domains of the nu- 
clear phase diagram. In heavy-ion (HI) collisions at in- 
termediate energies (bombarding energies between 20A 
MeV and 100,4 MeV) the projectile and target traverse 
each other in a time scale t < 50 fm/c forming a com- 
pressed and excited transient system that usually disas- 
sembles into several intermediate-mass-fragments (IMF) 
in a process known as "nuclear multifragmentation" pQ. 
The exact physical mechanism that drives this breakup 
is up to date an open issue. Two opposite scenarios have 
been considered: a fast breakup from an expanding source, 
consistent with a spontaneous spinodal mechanism, and a 
sequential slower breakup from a thermally equilibrated 
source. The knowledge of the relative importance of these 
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two mechanisms is essential to find-out a possible connec- 
tion between nuclear multifragmentation and the occur- 
rence of a "liquid-gas" -like phase transition [2] in nuclear 
systems with excitation energies in the range e*= (3 - 8)^4 
MeV, where a flattening of the "caloric curve" , signaling a 
possible liquid-gas phase transition, has been observed 
To discern between both possible mechanisms it is 
necessary to establish the breakup time of the multifrag- 
menting system. 

The experimental observable used in the present anal- 
ysis in order to probe the thermodynamical state of the 
produced transient nuclear systems, and in particular, to 
delimit the nuclear breakup time is thermal hard photon 
emission. At intermediate-energy HI collisions, photons 
above 30 MeV (hard photons) issue from bremsstrahlung 
in incoherent proton-neutron collisions jZj . The main hard 
photon contribution, the so-called "direct" component, is 
emitted from first-chance nucleon-nucleon (NN) collisions 
in the preequilibrium stage of the reaction. A second and 
softer flash of hard photons is emitted from secondary 
NN collisions in the later stages [HK| from a thermalizing 
source |10) . Hence, the experimental hard photon spec- 
trum measured in HI reactions has been found to be well 
reproduced by the sum of two exponential distributions 
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£L = K d e- E ^ E °+K t e- E ^ E ° (1) 

where d and t stand for direct and thermal, respectively, 
and the factors K^,t are related to the intensity of each 
source. The power of first chance bremsstrahlung photons 
is the fact that they result from incoherent collisions. This 
allowed one to equate the observed intensity and slope of 
the gamma spectrum with the number and hardness of 
collisions respectively. With the kinematic thermal model 
this information can also be obtained from secondary 
collisions; the soft spectrum now corresponds to the num- 
ber of collisions integrated over the lifetime of the equi- 
librated source and its formation, and the slope param- 
eter is related to the average temperature of the source 
[T2"] . Although multiple collisions lead to bremsstrahlung 
quenching (see e.g. J2|) ^ is believed not to play an impor- 
tant role because we exploit the photon spectrum above 
E lcu t=30 MeV. To observe about 10% quenching the col- 
lision rate of the nucleons in the thermal source should 
not exceed -E 7CU t/hbar=1.5 ■ 10 _1 (fm/c) _1 , in the model, 
without considering Pauli blocking, the rate typically is 
1.3 • 1CP 1 (fm/c) -1 . Thus the time scales derived in this 
work are not limited by the quenching phenomena. 

The analysis presented in this paper has two different 
goals. First, we want to elucidate the dependence of the 
thermal photon temperature on the reaction centrality in 
129 Xe + ™ Qt Sn reactions at 50,4 MeV; at this energy, this 
heavy and symmetric system provides good conditions for 
observing the effects of a possible nuclear liquid-gas phase 
transition, and moreover, the charged-particle distribu- 
tions and related observables have been studied in detail 
for this reaction by the INDRA Collaboration (see e.g.|14l 
I15U16II17II18| and references therein). The second objective 
consists in assessing the nuclear breakup time-scales us- 
ing the measured second-chance photon multiplicities and 
a thermal model that describes well the observed spectra. 
This investigation is performed for 7 different reactions 
studied by the TAPS Collaboration: 86 Kr + 58 Ni at 60.4 
MeV 0, 181 Ta + 197 Au at 40,4 MeV 0, 208 Pb + 197 Au 
at 30 A MeV jHj, 36 Ar + 107 Ag at 60 A MeV [H], 36 Ar + 
197 Au at 60,4 MeV [TO], 36 Ar + 58 Ni at 60,4 MeV HDJ and 
129 Xe + nat Sn at 50,4 MeV. The paper is organized as 
follows: in Sect. [3 the experimental setup is described. In 
Sect. |2| the particle identification methods are presented. 
The results obtained from the analysis of the inclusive 
hard photon spectrum are summarized in Sect.^J In Sect. 
[5] we focus on the Xe + na 'Sn exclusive measurements; 
first, we present the experimental centrality classes as well 
as the method employed to estimate the average impact 
parameter; next, we present the direct and thermal hard 
photon measurements for each centrality class. In Sect. El 
the thermodynamical properties extracted from the ther- 
mal hard photon signal are presented. In particular in 16. 21 
the results on the lifetimes of the thermal hard photon 
sources obtained for the different reactions are presented. 
Section gives a summary. 



2 Experimental setup 

The experiment E-300 was performed at the "Grand 
Accelerateur National d'lons Lourds" (GANIL) located in 
Caen, France, during 13 days of total beam time in May 
1998. The GANIL accelerator system delivered a Xe 46+ 
beam of 50,4 MeV with a bunch rate of 10 MHz and a 
nominal intensity of 5 nA, around 68 ions per pulse. The 
Xe ions impinged on a nat Sn target of 1 mg/cm 2 thick- 
ness (with associated interaction probability of 1.9T0 -3 
per burst). The experimental setup consisted of the "Two 
Arms Photon Spectrometer" TAPS associated with three 
particle multidetectors: the GANIL Silicon Strip Detector 
(SSD), the Washington University "Dwarf Ball" (DB) and 
the KVI "Forward Wall" (FW). A general layout of the 
experimental setup is presented in Fig. ^ 

The TAPS electromagnetic calorimeter ^1 consists of 
384 scintillation modules. Each one consisted of a BaF2 
hexagonal crystal, 2.95 cm inner radius and 25 cm length 
(corresponding to 12Ao radiation lengths) |20U21j . and a 
Charged Particle Veto (CPV) detector of plastic (NE102A) 
scintillator. TAPS offers a large solid angle coverage allow- 
ing for high quality angular distribution measurements in 
a wide energy range, from statistical (E 7 > 5 MeV) to 
the hardest photons (E 7 w 200 MeV). In this experiment 
TAPS modules were arranged in six blocks of 64 modules 
(8 x 8), which were placed in an almost symmetrical con- 
figuration in the horizontal plane around the target, and 
at an average distance of 56 cm, covering about 20% of 
the full solid angle. 

The GANIL silicon strip detector telescope SSD was 
placed at 19 cm downstream the target inside the vacuum 
chamber, covering the angular range of 2.2° < 8 < 10.3°. 
This device, consisting of a first disk 150 /im thick with 64 
semicircular strips and of a second disk 500 /im thick with 
128 radial strips, was sensitive to projectile-like fragments 
(PLF) emitted in binary reactions and intermediate-mass 
fragments (IMF) emitted in the forward direction from 
more dissipative reactions. 

The Washington University "Dwarf Ball" (DB) 
covered an angular range of 31° < 8 < 168° around 
the target. This detector system allows to identify light- 
charged particles (LCP) and IMF. The Dwarf Ball (DB) 
consists of 64 phoswich detectors forming a sphere with 
an inner radius of 41.5 mm. Each phoswich detector is 
made of a thin plastic scintillator and of a CsI(Tl) in- 
organic crystal scintillator. The plastic scintillator type 
(Bicron BC400 or BC446) depends on the covered angle 
and so does its thickness, 10 ^m and 40 fim for the most 
backward and forward modules, respectively. The CsI(Tl) 
crystal glued on the back of each plastic is 4 mm to 8 mm 
thick, depending also on the angle. 

The KVI "Forward WalF'(FW) |23], consisted of 92 
plastic scintillator phoswiches, detects and identifies LCP 
and IMF (up to Z«15) emitted in the forward direction. 
In our experimental setup, this detector was placed down- 
stream 74 cm away from the target covering the forward 
hemisphere behind the SSD, 2.5° < 8 < 25° and the 
whole azimuthal range. Each FW phoswich is composed 
of two organic plastic scintillators which are heat-pressed 
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TAPS blocks 




Fig. 1. General layout of the experimental setup used to study the reaction 129 Xe+ nat Sn at 50A MeV at GANIL. The detector- 
system consists of the Silicon Strip Detector (SSD), the Washington University Dwarf Ball (DB), the TAPS electromagnetic 
calorimeter and the KVI Forward Wall (FW). 



together; a 1 mm thick "fast" (r = 2.4 ns) NE102A plas- 
tic followed by a 50 mm thick "slow" (r = 320 ns) NE115 
plastic. 

This complete experimental setup allowed correlated 
measurements of photons and charged-particles, necessary 
to study the dependence of the hard-photon production on 
the reaction topology. The photon data were recorded with 
the minimum bias trigger for neutral particles (7*DB), 
which signaled events with at least one neutral hit in 
TAPS of energy > 10 MeV detected in coincidence 
with one or more charged particles in the Dwarf Ball. The 
minimum bias (MB) reaction trigger was defined by the 
condition of detecting at least a charged particle in the 
Dwarf Ball. 



3 Particle identification 

Photons are identified in TAPS by means of a Pulse Shape 
(PSA) vs. time-of-flight (TOF) analysis and the CPV in- 
formation. Photons and leptons induce a stronger inten- 
sity of the fast BaF 2 light component compared to hadronic 
particles, producing therefore a higher ratio of the fast 
to the slow energy components in the BaF2 (PSA). Pho- 
tons and relativistic electrons are located in time within 
a prompt peak centered at TOF=1.87 ns (with a reso- 
lution a = 340 ps [U); these particles produce hence a 
recognizable contour in the PSA-TOF spectrum. In order 



to separate photons from relativistic electrons the infor- 
mation (fired or not) delivered by the CPV of each mod- 
ule is used. Nevertheless, given that cosmic muons have 
mainly vertical trajectories not firing the CPV, and have 
random TOF, they can be misidentified as photons. This 
is avoided by defining a PSA-TOF cosmic contour several 
nanoseconds far from the photon contour but with the 
same dimensions. Subsequently, the energy and angular 
spectra measured within this contour (see Fig. |2J) are sub- 
tracted from the raw photon energy and angular spectra 
respectively. 

The Z\E-E telescope technique [241125} allows the charge 
identification of LCP and IMF detected in the DB and in 
the FW, as well as of the PLF and IMF impinging on the 
SSD. The light measured in the first stage of each multi- 
detector (the first component of each phoswich in the DB 
and the FW, or the first disk in the SSD) and the sig- 
nal of the second component are plotted in bidimcnsional 
histograms in which the different curves correspond to dif- 
ferent Z. In the DB, additional LCP isotopic identification 
was also possible exploiting the pulse shape properties of 
the CsI(Tl) light. 



4 Inclusive hard photon spectra in 129 Xe + 
na *Sn at 50A MeV reactions 

Figure |21 shows the inclusive raw photon energy spectrum 
measured in the 129 Xe + ™ at Sn reaction from the collection 
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E y (MeV) 



Table 1. Summary of the experimental inclusive hard photon 
(E 1 > 30 MeV) results obtained for the 129 Xe + nat Sn reac- 
tion at 50 A MeV. For the direct and thermal 7 components 
we quote: inverse slope Eq, relative intensity I, source velocity 
(3, cross-section <r 7 , multiplicity M lt and (only for the direct 
component) hard photon probability P 7 . 



Direct 



Thermal 



So (MeV) 


15.6 ± 1.0 


7.0 ± 0.6 


/(%) 


78 ± 1 


22 ± 1 


P 


0.15 ± 0.01 


0.16 ± 0.01 


cr 7 (mb) 


4.9 ± 0.6 


1.4 ± 0.2 


M 1 


(9.3 ± 0.8)- 10" 4 


(2.6 ± 0.3)- 10" 4 


R f 


(7.8 ± 0.7) ■ 10~ 5 





Fig. 2. Inclusive raw photon spectrum measured in the center- 
of-mass frame for the reaction 129 Xe + nat Sn at 50 A MeV. 
The triangles show the cosmic ray background contribution. 



of 5.8 • 10 6 photons with E 1 > 10 MeV detected in TAPS 
and transformed into the CM frame ((3aa ~ Pnn ~ 0.15). 
After cosmic-ray subtraction, and for energies E 1 > 30 
MeV this spectrum is well described by the double source 
fit of Eq.^(Fig. 3). The direct component exhibits an in- 
verse slope parameter Eq = (15.6 ± 1.0) MeV whereas the 
thermal slope parameter is Eq = (7.0 ± 0.6) MeV. Both 
slopes are in agreement with the direct [26| and thermal 
12 photon slope systematics measured in intermediate 
A + A reactions. Thermal hard photons account for 22% 
of the total yield above 30 MeV. The strength and in- 
tensity of the thermal emission are also confirmed by a 
double source analysis of the hard photon angular distri- 
bution; however, due to the symmetry of the system, the 
NN and A A source velocities cannot be disentangled |27) . 

The direct (thermal) hard photon multiplicity, i.e. the 
number of direct (thermal) hard photons emitted per nu- 
clear reaction, has been obtained via the following expres- 
sion: 



Ml 



7 




E 
n 

£ 10" 



30 35 40 45 50 55 60 
E,(MeV) 



30 35 40 45 50 
E,(MeV) 



Fig. 3. Left panel: experimental inclusive hard photon spectrum 
measured for the reaction 129 Xe + nat Sn at 50A in the range 
i? 7 =30 - 65 MeV, fitted to the Eq. QJ the solid (dashed) line is 
the direct (thermal) exponential distribution. Right panel: in- 
clusive thermal hard photon spectrum obtained after subtracting 
the fitted direct component. 



where i stands for direct or thermal, ct 7 is the experimen- 
tal hard photon cross section determined by integrating 
the corresponding da fdE 1 spectrum above E 1 — 30 MeV, 
and <jr = (5300 ± 600) mb is the experimental total reac- 
tion cross section for Xe+Sn, obtained from the charged 
particle distribution measured in the DB with the MB 
reaction trigger. We report the direct and thermal hard 
photon cross sections and multiplicities in Table ^ The 
direct hard photon probability Pf?, i.e. the probability to 
produce a hard photon in a first-chance proton-neutron 
collision, is determined from the experimental direct hard 
photon multiplicity: 



pd 

7 



AE' 



<JR ■ (N pn ) b (N pn ) b 



(3) 



where (A p „)ti is the number of first chance pn collisions 
averaged over impact parameter, calculated from the geo- 
metrical "equal-participant" model of Nifenecker and Bon- 
dorf [2Hj- (N pn ) b is = 11.79 for the Xe on Sn system. The 
measured in 129 Xe + not Sn reactions at 50,4 MeV is 
(7.8 ± 0.7) ■ 10 -5 . We note that we do not quote a ther- 
mal hard photon probability, since we should determine 
the average number of secondary (and not primary) pn 
collisions, from which the thermal component originates, 
and this number is not well determined experimentally nor 
theoretically. More details of the inclusive analysis can be 
found in |2J. 
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Table 2. Conditions on the fragment multiplicity in the SSD (M SSD ), the charged particle multiplicity (IMF + LCP) in the 
DB (M ) and the FW (M FW ) that define each centrality class. The contribution to the reaction cross-section, and to the 
nuclear reactions detected by the minimum bias photon trigger ( r y*DB) for each centrality class are reported. 



Centrality class 


Multiplicity condition 




%o- jR ( 7 *DB) 


A or Periph. 


M DB+FW = x _ 2 M SSD { g <5)=0 


57% 


23% 


B 


M DB+FW = 12 


67% 


26% 


C 


2 < M DB+FW < 8 


38% 


58% 


D 


5 < M DB+FW < 9 


10% 


24% 


E 


6 < M DB+FW < 10 


7% 


18% 


F or Central 


9 < M DB+FW < 15 


2% 


6% 



Table 3. Characteristics of the six centrality classes considered for Xe+Sn reactions at 50A MeV. For each centrality bin we 
quote: average charged particle multiplicity M cp and impact parameter {b) geo estimated with the geometrical model of Cavata et 
al. \Hyf: direct photon multiplicity (M$), average number of first chance proton-neutron collisions (N pn ) obtained as described 
in ]3(M3tf . and associated averaged impact parameter (6) 7 estimated using the "equal-participant" model 



Centrality class 


(M ap ) 


(b) 9eo (fm) 








(N pn ) 


(b) 7 (fm) 


A or Periph. 


1.5 


9.4 ± 0.8 


(3.2 ± 0.2) 


10" 


-4 


4.1 ± 0.6 


9.2 ± 0.3 


B 


1.5 


9.4 ± 0.8 


(4.1 ± 0.3) 


10" 


-4 


5.2 ± 0.6 


8.8 ± 0.3 


C 


4.7 


5.2 ± 0.4 


(1.5 ± 0.1) 


10" 


-3 


18.6 ± 2.2 


5.9 ± 0.5 


D 


6.8 


3.3 ± 0.3 


(2.4 ± 0.2) 


10" 


-3 


30.7 ± 3.7 


4.0 ± 0.6 


E 


7.8 


2.6 ± 0.2 


(2.8 ± 0.2) 


10" 


-3 


36.1 ± 4.4 


3.1 ± 0.7 


F or Central 


11.4 


0.9 ± 0.1 


(3.8 ± 0.4) 


10" 


-3 


48.9 ± 6.7 


1.3 ± 1.1 



5 Exclusive hard photon analysis in 129 Xe + 
nat Sn at 50,4 MeV reactions 

5.1 Selection of centrality classes 

Six centrality classes of 129 Xe + raat Sn collisions have been 
selected to investigate the dependence of the thermal hard 
photon production on the reaction centrality. This choice 
results from a compromise between requiring a wide im- 
pact parameter coverage and obtaining enough hard pho- 
ton statistics for each centrality class. We report in Table 
13 the conditions on the fragment multiplicities measured 
in the SSD, Dwarf Ball and Forward Wall that define the 
selected centrality bins. 



measured and from the P^, which is considered to de- 
pend only on the Coulomb corrected bombarding energy. 
The obtained (N pn ) values are related to the impact pa- 
rameter b by the geometrical "equal-participant" model of 
of Nifenecker and Bondorf [2Hj • We note that if the impact 
parameter is obtained directly from the total hard photon 
yield, as commonly done, instead of only from the direct 
component, the deduced impact parameter might be dis- 
torted for reactions where a thermal photon contribution 
cannot be neglected. The average impact parameter for 
each reaction class, obtained by means of the direct hard 
photon multiplicities and by means of the charged particle 
multiplicities measured in the DB, are reported in Table 
El We find a good agreement between both methods, and 
in addition the whole impact parameter range is covered. 



5.2 Estimation of the impact parameter for minimum 
bias data 

In order to estimate the impact parameter b in 129 Xe + 
raat Sn all reactions, recorded by means of the MB trigger, 
we apply the geometrical method proposed by Cavata et 
al. j2Hl, which relies on the monotonous decrease of Mcp 
as a function of b. An alternative method to estimate the 
impact parameter in experiments where hard photon emis- 
sion can be analyzed is based on the correlation of the 
direct hard 7 yield with the impact parameter |30II31| . In 
this method, for a given reaction class, the experimental 
average number of first chance pn collisions (N pn ) is es- 
timated, making use of Eq. from the experimentally 



5.3 Determination of the impact parameter for 
reactions where a hard photon is produced 

We want to exploit the properties of the thermal hard pho- 
ton component as a tool to extract the thermodynamical 
state of the radiating source. It is therefore of crucial im- 
portance to determine as accurately as possible any bias 
introduced in the estimation of the impact parameter for 
particular reactions where a hard photon is produced. 

Within a given centrality class, defined by a range in 
the M cp distribution, hard photons are mostly emitted in 
reactions with high M cp values, characterized by a large 
projectile-target overlap, and consequently by a high num- 
ber of pn collisions. This bias through lower impact pa- 
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rameters can be reproduced by the geometrical "equal- 
participant" model; the reaction cross-section as a func- 
tion of b for reactions where a hard photon is produced, 
cr^(b), is modulated by the product of crfi(b) and N pn (b). 
As a consequence, the ratio between the reaction cross- 
section crjj(b) and the total (integrated over b) er^ is a 
representable function sensitive to the cross section devi- 
ation through lower b (see Fig. 0}. It is thus possible to 
estimate the average impact parameter for reactions where 
a hard photon is produced by means of an extension of the 
geometrical model of Cavata et al., correlating the exper- 
imental a^Mcp) I ijft and the theoretical N pn ) / 'ctr 
obtained from the geometrical "equal-participant" model. 
The average b of each centrality class has been obtained 
with this method after weighting the b(M cp ) by the con- 
tribution of each M cp within the centrality class. Table 0] 
reports the (b) estimated for each RC for reactions where 
a hard 7 is emitted. 



Table 4. Estimated average impact parameter (b) for each of 
the six centrality classes considered. 



Centrality class 


(b) 9eo (fm) 


Class A or Peripheral 


7.6 ± 0.3 


Class B 


7.2 ± 0.3 


Class C 


4.2 ± 0.4 


Class D 


2.8 ± 0.4 


Class E 


2.2 ± 0.6 


Class F or Central 


1.0 ± 0.9 




b(fm) 



Fig. 4. Evolution of on(b)/oR as a function of the impact pa- 
rameter without (empty histogram) and with (filled histogram) 
the condition that a hard photon is produced. 



5.4 Centrality dependence of the hard photon 
production 

The two-exponential fit of Eq.^ reproduces the hard pho- 
ton energy spectrum of each one of the six centrality classes, 
from peripheral to central events. The corresponding di- 
rect and thermal slopes parameters and intensities are 
listed in Table |5| Whereas the measured E^b) differ by 
up to 24% from the inclusive value, the direct slope vari- 
ation amounts to less than 5% of the inclusive Eq. The 
measurement of an almost constant direct slope confirms 
that the E$ is an observable that depends just on the 
bombarding energy [U] . Only for very peripheral reactions 
it is expected to observe a decrease of the direct photon 
slope mainly due to the decrease of the Fermi momentum 
in the surface of the nucleus [3T] . In Fig. the values of 
the thermal (direct) slopes are displayed as a function of 
the thermal (direct) multiplicities. The direct hard pho- 
ton slope is independent of the direct hard photon mul- 
tiplicity. This shows that the increase of the direct hard 
7 multiplicity with the centrality is just a consequence of 
the associated increase of the number of first-chance pn 
collisions. On the other hand, the thermal slope exhibits a 
linear dependence on the thermal multiplicity. This trend, 
observed for the first time, indicates that, assuming that 
the size of the thermalized source remains almost con- 
stant with centrality, the thermal production is sensitive 
to the excitation of the nuclear system. This result is sup- 
ported by the inclusive thermal systematics collected by 
the TAPS Collaboration: the thermal slopes and thermal 
multiplicities divided by the size of the system, scale with 
the energy available in the nucleus-nucleus center-of-mass 

5.5 Exclusive analysis of the GDR emission 

Exploiting the wide dynamic range of TAPS, we have 
estimated the 7 yield from the Giant Dipole Resonance 
(GDR) decay emitted in peripheral 129 Xe + na 'Sn reac- 
tions. The main motivation of this investigation was to 
identify a possible thermal bremsstrahlung and GDR mix- 
ture in the low energy range (E 7 ss 30 McV) of the hard 
photon spectra, which could distort the measurement of 
the thermal hard photon yield. Peripheral reactions ex- 
hibit the highest probability of GDR7 contamination be- 
cause the reached excitation energies can still be not high 
enough for suppression of the collective mechanisms re- 
sponsible of the GDR7 emission > and because of the 
softer observed thermal emission. So the GDR analysis in 
peripheral reactions allowed us to conclude about a pos- 
sible GDR7 contamination for all reaction classes. 

In order to perform the GDR7 analysis, the direct and 
thermal hard photon components have been subtracted 
from the total photon spectra. This subtraction is done 
by extrapolating the double hard 7 exponential fit down 
to the lowest energy, given by the TAPS LED threshold of 
10 MeV. This extrapolation technique is routinely applied 
in GDR measurements from HI reactions (see e.g. j-i2U3.il 

ED). 



R. Ortega et al.: Constraints on the time-scale of nuclear breakup from thermal hard-photon emission 



7 



Table 5. Characteristics of the hard photon spectra measured in the six different centrality classes: direct (Eq ) and thermal 
(Eq) slopes, ratio of thermal to total hard-photon intensities (h/ltot), and thermal photon multiplicity (Ml). 



Centrality class E% (MeV) E* (MeV) I t /Itot(%) M* 



X 2 /ndf 



A or 


Periph. 


14.9 ± 0.9 


5.7 ± 0.5 


19 ± 1 


(7.6 ± 0.7)40" 


5 


1.5 


B 




14.6 ± 0.9 


5.7 ± 0.5 


19 ± 1 


(8.9 ± 1.1)-10~ 


5 


1.2 


C 




15.6 ± 0.9 


6.3 ± 0.6 


19 ± 1 


(3.4 ± 0.4)-10~ 
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Fig. 5. Direct and thermal hard photon slopes vs their respec- 
tive multiplicities for the six Xe+Sn centrality classes consid- 
ered. 



We have compared the GDR spectra measured for pe- 
ripheral centrality classes A and B with the GDR dis- 
tributions measured with the multidetector MEDEA and 
reported by T. Suomijarvi et al. |34| . In this work, the 
GDR7 emission from hot nuclei of A 115 formed in the 
36 Ar + 90 Zr at 27,4 MeV reaction is investigated in de- 
tail. The comparison of the GDR7 yield measured in this 
analysis and the one measured for peripheral reactions in 
our work comes easily to mind when considering the sim- 
ilarities in mass and in excitation energy of the emitting 
sources: 

— In peripheral 129 Xe + nat Sn reactions the quasi-projectile 
and quasi-target are close to the 129 Xe and na 'Sn masses. 
Assuming that both, quasi-target and quasi-projectile, 
develop collective GDR oscillations after the collision, 
the photon spectrum measured in the peripheral ccn- 
tralities scaled by a factor of ~l/2 can be compared 
to the GDR 7 yield emitted by a nuclei of A ~ 115; 



In spite of the different geometry and bombarding en- 
ergy of both reactions, the excitation energies reached 
at large impact parameters in 129 Xe + ™ at Sn reactions 
(~ 3,4 MeV, see next section) are comparable with the 
e* attained in the 36 Ar + 90 Zr MeV reaction at 27 A. 
Besides, the GDR yield measured for the 36 Ar + 90 Zr 
system has been found to remain constant with the 
excitation energy jM] (see Fig. 0. 
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Fig. 6. Spectra in the GDR region for the 36 Ar + 90 Zr reaction 
\34\ , for two different excitation energies: e* =350 MeV (dots), 
e* = 500 MeV (squares) ]?!^ ! compared to the experimental 7 
spectrum obtained in peripheral 129 Xe + nat Sn reactions (tri- 
angles) after subtraction of the hard-photon yield (extrapolated 
down to 10 MeV). 



The GDR yields (the 7 spectra after subtraction of the 
bremsstrahlung contribution) for the 36 Ar + 90 Zr reaction 
for e* = 350 MeV and for e* = 500 MeV from gl, and the 
GDR yield estimated for the 129 Xe + nat Sn centrality class 
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A 1 are plotted together in Fig. El The global compatibility 
of the three distributions assures the validity of the esti- 
mation of the GDR yield for the peripheral Xe+Sn hard 
photon spectra. The difference found in the 25-30 MeV 
region between the TAPS and MEDEA measurements is 
due to the fact that in the latter the subtracted hard pho- 
ton distribution is parameterized by a single exponential. 
Taking into account that for peripheral reactions the to- 
tal photon multiplicity is around 9 -10 -5 at E 1 — 30 MeV 
and that the GDR yield at this energy is found to be lower 
than 10~ 5 (see Fig. |SJ, we conclude that the GDR con- 
tamination in the low energy range of the hard photon 
spectra can be neglected. This result is in good agreement 
with the highest energy GDR gamma contribution pre- 
dicted from the systematics. Indeed, taking the centroid 
energy of the GDR as Egdr = 76. 5A -1 / 3 and the max- 
imum width r GDR w 12 MeV |53I55] . the GDR emission 
neither from the Xe projectile nor from the Sn target are 
expected to exceed 28 MeV. 



6 Derived (thermo)dynamical properties 



Table 6. Nuclear temperatures T estimated through Eq.^from 
the experimental thermal hard photon slopes Eq of the inclusive 
and exclusive hard photon spectra measured in 129 Xe + nat Sn 
at 50 A MeV. 



Centrality class 


El (MeV) 


T (MeV) 


Inclusive 


7.0 ± 0.6 


5.5 ± 0.8 


A or Peripheral 


5.7 ± 0.5 


4.4 ± 0.4 


B 


5.7 ± 0.5 


4.4 ± 0.4 


C 


6.3 ± 0.6 


4.9 ± 0.5 


D 


7.8 ± 0.7 


6.1 ± 0.5 


E 


7.9 ± 0.7 


6.2 ± 0.5 


F or Central 


8.7 ± 0.8 


6.8 ± 0.6 



6.2 Lifetime of the thermalized system 

The multiplicity of hard photons emitted by an equili- 
brated nuclear fragment can be calculated by means of 
the model of Neuhauser and Koonin from the following 
expression [T2] : 



6.1 Nuclear Temperature 

Although, due to their long mean free path in nuclear mat- 
ter, thermal hard photons are not in thermal equilibrium 
with the source, their spectral shape is directly correlated 
to the temperature of the radiating nuclear fragment. As 
in our previous analyses |12| . we use the kinetic thermal 
model of Neuhauser and Koonin JT] to calculate the re- 
lation between the inverse slope Eq and the temperature 
T of the nuclear system. The electromagnetic radiation 
rate emitted by a hot and equilibrated nuclear fragment 
is expressed in this model as a function of the local den- 
sity and temperature of the source. It is assumed that 
bremsstrahlung in proton-neutron collisions is the main 
source of thermal photon emission for E~ > 30 MeV. The 
hard photon spectra thus calculated are accurately ap- 
proximated by an exponential with slope Eq in agreement 
with the experimental data ^2]. The evolution of Eq val- 
ues extracted from a fit above 30 MeV with the temper- 
ature, T, is quantitatively well described by the following 
linear expression: 



T(MeV) = (0.78 ± 0.02) • E^MeV) 



(4) 



in the region T w 3-10 MeV and p « (0.3-1.2)p P] , 
We have applied this photon thermometer to extract the 
temperature of the radiating system produced in the dif- 
ferent 129 Xe + ™ at Sn centralities studied in this work. The 
obtained temperatures are reported in Table El The tem- 
perature of the nuclear system for central collisions (cen- 
trality class F, b/b max < 0.1) is a 50% ± 20% higher than 
for peripheral collisions (centrality class A, b/b max w 0.6). 



the GDR yield measured for centrality classes A and B are 
almost identical, so only the yield for class A is shown in the 
figure. 



dt 



.in 



dR» K (T,p) jri 



(5) 



where R^ K (T,p) is the rate of photons emitted from a 
thermal nuclear fragment with temperature T and den- 
sity p. From this model, R^ K (T,p) is found to scale with 
~ T 6 ' 7 and to be proportional to p. The quadratic depen- 
dence of the photon rate on p expected naively turns out 
to be linear due to the Pauli blocking. The photon rate 
can be hence expressed as R^ K ~ Rq K - T e T ■ p, with the 
constant R VK = 1.1 • 10" 13 MeV^ ^fm/cT 1 . We can 
further simplify Eq. [SJ by assuming that the temperature 
is uniform in the volume V and roughly constant during 
the emission time, 



V ■ At-R, 



NK 



T 



6.7 



' P 



(6) 



where V and At are respectively the volume and the life- 
time of the radiating nuclear source. The volume can be 
approximated as the ratio of the sum and the projectile 
and target nucleons 2 over the nuclear density, 
V«(A t + A p )/p. 

With these simplifying assumptions, the thermal hard 
photon yield, scaled to the relative size of the system, can 
be written as: 



M, 



NK 



(A t + A p ) 



At-R, 



NK 



T 



6.7 



(7) 



This relation, which is independent of the density, can be 
exploited to estimate the lifetime At of the thermal equi- 
librated source, by fitting Eq. (7) to the experimentally 
measured multiplicities (reported in Table 5) as a function 

2 We neglect here emission of particles during preequilib- 
rium. 
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of the temperature. Figure shows the experimental (in- 
clusive) Mi, divided by size of the system A to t — A t +A p , 
as a function of the estimated T, for the Xe+Sn and the 
six other reactions studied |T2]. The lines correspond to 
the calculated M^ K as a function of the T for two differ- 
ent lifetimes (Z\r=100 fm/c and Z\r=35 fm/c). 

The calculated multiplicities successfully reproduce the 
experimental thermal hard photon yields measured for all 
the reactions (see Figs. and 0, with At values about 
100 fm/c, which are consistent with the expected lifetime 
of an equilibrated source (see e.g. [HE])- This result con- 
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Fig. 7. Experimental thermal bremsstrahlung multiplicity, M~, 
divided by the size of the colliding system, Atot ~At+A p , plot- 
ted as a function of the extracted nuclear temperatures, for 
the different reactions studied by the TAPS Collaboration. The 
lines correspond to two different values of the source lifetime 
estimated with Eq. ^ 



firms the validity of the thermal model, since it reproduces 
both experimental observables, the slope of the thermal 
hard photon spectrum and the thermal multiplicity, with 
the temperature as the single input parameter. In Figs. 
and 03 it can be observed that all the systems studied 
have a lifetime At of the order or larger than 100 fm/c, 
except the two lighter systems at 60,4 MeV, 86 Kr+ 58 Ni 
and 36 Ar+ 58 Ni, which seem to survive a shorter At ~ 35 
fm/c, similar to the transit time of the colliding ions. In 
these two reactions, the achieved excitation energy is also 
higher than in the rest of reactions fT^] . This difference in 
At might indicate that equilibrated nuclear sources break 
faster when they have a smaller size or higher excitation 
energies, as observed in [5rjll37U38| . The experimental ther- 
mal hard photon multiplicities measured for the 129 Xe + 
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Fig. 8. Thermal hard photon multiplicities, Ml, divided by the 
total size of the system, Atot =A t +A p , plotted as a function of 
the nuclear temperature extracted from Eq. fj] for the different 
129 Xe + nai Sn centrality classes. The lines correspond to two 
different values of source lifetime computed with Eq. ^ 



na 'Sn centrality classes exhibit a trend compatible with 
lifetimes of the order of 100-300 fm/c (FigJ^J. In our mea- 
surements, within the data uncertainty, we do not observe 
any change of the lifetime of the equilibrated fragment 
with the impact parameter of the collision for Xe+Sn. 
The relatively large lifetime values obtained, jointly with 
the measurement of an important production of thermal 
hard photons, suggest that, at least for the investigated 
Xe+Sn reactions, the fireball does not undergo an instan- 
taneous breakup from a rapidly expanding state (spinodal 
decomposition). 



6.3 Caloric curve constructed from Xe + Sn exclusive 
measurements 

As done in a previous paper for inclusive thermal pho- 
ton measurements |12) . we have measured the relation- 
ship between the estimated temperature and the excita- 
tion energy, e*, of the hot equilibrated fragment formed 
in Xe + Sn reactions at the different centralities studied 
in this work. At variance with experiments that measure 
the excitation energy with charged particle detectors, we 
cannot determine "calorimetrically" e* from our data. The 
excitation energy of Xe quasiprojectiles produced from pe- 
ripheral to central collisions has been deduced by Steck- 
meyer et al. |18 | as a function of the charged particle mul- 
tiplicity measured by the INDRA detector. The methods 
described in [TBI are based on the determination of the 
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velocity of the quasiprojectile. The excitation energy is 
calculated from the kinetic energies of products belonging 
to the source. We have employed the geometric method 
of Cavata et al. |29j to relate the average impact parame- 
ter of each Xe + Sn reaction class with the corresponding 
value of M ' C p DRA 1 in order to obtain the measured value of 
e* for each centrality class. The estimated M™ DaA and e* 
values are reported in Table \7\ 



Table 7. Values of the charged particle multiplicity measured 
with INDRA, Mcp DnA , and of the estimated excitation energy, 
e* , for each Xe+Sn centrality class with averaged impact pa- 
rameter (b). 



> 



10 



8- 



4- 



2 t/ 



+ 



Centrality class 


(b) 


(fm) 


M' C » DRA 


e*(A MeV) 


Inclusive 


3.8 


± 


1.0 


25 ± 1 


5.6 ± 1.0 


A or Peripheral 


7.6 


± 


0.3 


14 ±1 


2.8 ± 0.3 


B 


7.2 


± 


0.3 


16 ± 1 


3.3 ± 0.4 


C 


4.2 


± 


0.4 


24 ± 1 


5.2 ± 0.6 


D 


2.8 


± 


0.4 


27 ± 1 


6.1 ± 0.7 


E 


2.2 


± 


0.6 


29 t\ 


6.7 ± 0.7 


F or Central 


1.0 


± 


0.9 


32 +_% 


7.7 t\i(8.7 ± 1.5) 



10 12 
e*(A MeV) 



Fig. 9. Caloric curve constructed from the thermal hard pho- 
ton slope measurements for the 129 Xe + nat Sn reaction at 50 A 
MeV (Table 6) and the estimated excitation energies (Table 7), 
compared to the Fermi liquid model curves with inverse level- 
density parameters K—8/A MeV (dashed line) and K=13/A 
MeV (dotted line). 



The errors of e* take into account the uncertainty of (b) 
and a ~10% systematic error due to the reconstruction of 
the quasiprojectile (T^l- As expected, the excitation energy 
decreases linearly with increasing impact parameter. Since 
the highest measured charged particle multiplicity does 
not correspond to a real (b) = 0, the geometrical method 
used to relate (b) and M™ DRA is not accurate enough for 
central collisions. We therefore also quote in parentheses 
for centrality class F the value of e* predicted for central 
collisions by the model of Natowitz et al. [HHI • The caloric 
curve obtained is displayed in Fig. [5] The temperature 
increases smoothly with the excitation energy, the (e*,T) 
pairs o btain ed are below the expected Fermi fluid curves 
(T = /Re* (MeV) with K = 8/ A - 13/ A MeV covering 
the known range of average nuclear level density parame- 
ters 00]). Peripheral reactions (classes A and B) lead to 
excitation energies e*« 3^4 MeV, near those identified as 
the onset of the leveling region of the caloric curve for 
nuclear systems of size A= 180-241 0]. 

7 Summary 

Hard photon [E~ > 30 MeV) emission has been stud- 
ied in 129 Xe + ,mt Sn at 50,4 MeV for six bins of impact- 
parameter. The thermal component has been determined 
from a double source analysis of the measured spectra. 
For each centrality class, the temperature of the produced 
nuclear systems has been extracted from the measured 
thermal hard photon slope, by applying a thermal nu- 
clear bremsstrahlung model that successfully reproduces 
the measured thermal hard photon spectral shape and 
multiplicities. The temperature obtained exhibits a small 



but systematic increase with decreasing impact param- 
eter, from T — 4.4 MeV for the most peripheral colli- 
sions (b/bmax ~ 0.6) to T = 6.8 MeV for the most cen- 
tral class ( b/bmax < 0.1). The amount of second-chance 
bremsstrahlung emission measured in each centrality bin 
is well described by the thermal model by assuming emis- 
sion times of the order r 100 - 300 fm/c for all central- 
izes. The lifetimes of the produced equilibrated sources 
have also been estimated for six other heavy-ion reactions 
for which thermal bremsstrahlung emission has been mea- 
sured by the TAPS Collaboration. The inclusive thermal 
hard photon multiplicities measured as a function of the 
temperature T for the different reactions suggest lifetimes 
r ~ 100 fm/c, i.e. between two and three times longer 
than the transit time of the colliding ions at these ener- 
gies. Such relatively large time scale disfavor scenarios of 
instantaneous spinodal break-up and is more consistent 
with sequential fragmentation from a thermally equili- 
brated source. It is worth noting, however, that for the two 
smallest target-projectile combinations, 86 Kr+ 58 Ni and 
36 Ar+ 58 Ni at 60A MeV, which are the systems with more 
energy available in the nucleus-nucleus center-of-mass, the 
experimental thermal photon multiplicities seem to be bet- 
ter reproduced with shorter lived equilibrated sources (r ~ 
35 fm/c). 

Finally, we present for the first time a caloric curve, 
e*(T), obtained from thermal hard photon measurements 
for different impact parameter bins of the same reaction, 
Xe -(- Sn. Such a curve falls below the one expected for 
a Fermi liquid, in agreement with the curve constructed 
from thermal hard photon slopes measured in different 
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heavy-ion reactions ^2], and consistent also with the avail- 
able systematics of nuclear caloric curves 



8 Acknowledgements 

We thank the GANIL accelerator staff for providing a high 
quality beam and technical support during the experi- 
ment. We thank J.C. Steckmeyer and the INDRA Col- 
laboration for providing the charged-particle multiplicity 
data, and the MEDEA Collaboration for providing the 
GDR data. This work has been in part supported by 
IN2P3-CICYT PN97-1 agreement, by the Dutch Foun- 
dation FOM, by the European Union HCM network un- 
der Contract No. HRXCT94066, and by DGICYT under 
contract FPA2000-2041-C02-01. R. Ortega aknowledges 
support given by the Alexander von Humboldt Founda- 
tion. D.d'E. acknowledges support by the European Union 
TMR Programme (Marie-Curie Fellowship No. HPMF- 
CT-1999-00311). 



29. C. Cavata et al, Phys. Rev. C294, 47 (1993). 

30. S. Riess et al, Phys. Rev. Lett. 69, 1504 (1992). 

31. G. Martinez et al, Phys. Lett. B334, 23 (1994). 

32. J.J. Gaardhoje et al, Phys. Rev. Lett. 59, 1409 (1987). 

33. J.H. Le Faou et al, Phys. Rev. Lett. 74, 3321 (1994). 

34. T. Suomijarvi et al, Phys. Rev. C53, 2258 (1996). 

35. B.L. Berman and S.C. Fultz, Rev. Mod. Phys. 47, 713 
(1975). 

36. R. Wada et al, Phys. Rev. C69, 044610 (2004). 

37. L. Beaulieu et al, Phys. Rev. C84, 5974 (2000). 

38. Z. He et al, Phys. Rev.C63, 011601 (2001). 

39. J.B. Natowitz et al, Z. Phys. A325, 467 (1986). 

40. B. Nerlo-Pomorska et al, Int. J. Mod. Phys. E 14 (2005) 
505 



References 

1. J. P. Bondorf, A.S. Botvina, A.S. Ilinov, I.N. Mishustin and 
K. Sneppen, Phys. Rept. 257 (1995) 133; L.G. Moretto and 
G.J. Wozniak, Ann. Rev. Nucl. Part. Sci. 43 (1993) 379. 

2. V. E. Viola, Nucl. Phys. A 734 (2004) 487. 

3. J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995). 

4. J. B. Natowitz, Phys. Rev. C65, 034618 (2002). 

5. J. B. Natowitz et al, Phys. Rev. L ett. 89, 212701 (2002). 

6. W. Trautmann, |micl-ex/0411023| Proceedings INPC2004, 
Goeteborg, Sweden, June 27- July 2 (2004). 

7. H. Nifenecker and J. Pinston, Ann. Rev. Nucl. Part. Sci. 
40, 113 (1990). 

8. G. Martinez et al, Phys. Lett. B349, 23 (1995). 

9. Y. Schutz et al, Nucl. Phys. A622, 404 (1997). 

10. D. d'Enterria et al, Phys. Rev. Lett. 87, 22701 (2001). 

11. D. Neuhauser and S.E. Koonin, Nucl. Phys. A462, 163 
(1987). 

12. D. d'Enterria et al, Phys. Lett. B538, 27 (2002). 

13. M.J. van Goethem et al, Phys. Rev. Lett. 88, 122302 
(2002). 

14. N. Marie Phys. Lett. B391, 15 (1997). 

15. N. Marie et al, Phys. Rev. C58, 256 (1999). 

16. D. Gourio et al, Eur. Phys. J. A7, 245 (2000). 

17. N. Bellaize et al, Nucl. Phys. A709, 367 (2002). 

18. J.C. Steckmeyer et al, Nucl. Phys. A686, 537 (2001). 

19. R. Novotny, IEEE Trans. Nucl. Sci. 38, 379 (1991). 

20. M. Laval et al, Nucl. Inst, and Meth. A206, 69 (1983). 

21. A. Gabler et al, Nucl. Inst, and Meth. A346, 164 (1994). 

22. D.W. Stracener et al, Nucl. Inst, and Meth. 294, 485 
(1990). 

23. S. Luke et al, Phys. Rev. C47, 1211 (1993). 

24. R. Ortega, Ph.D.Thesis, Univ. Aut. de Barcelona, 2003. 

25. L. Tassan-Got, Nucl. Inst. Meth. B194, 503 (2002). 

26. J. van Pol, Ph.D.Thesis, Rijkuniversiteit Groningen, 1995. 

27. R. Ortega, F. Fernandez, D. D'Enterria and G. Martinez, 
Acta Phys. Polon. B 33, 965 (2002). 

28. H. Nifenecker and J. P. Bondorf Nucl. Phys. A442, 478 
(1985). 



